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ABSTRACT 

A long standing prediction of steady state gas-phase chemical theory is that 
H2O and O2 are important reservoirs of elemental oxygen and major coolants of the 
interstellar medium. Analysis of SWAS observations has set sensitive upper limits 
on the abundance of O2 and has provided H2O abundances toward a variety of star 
forming regions. Based on these results, we show that gaseous H2O and O2 arc not 
dominant carriers of elemental oxygen in molecular clouds. Instead the available 
oxygen is presumably frozen on dust grains in the form of molecular ices, with a 
significant portion potentially remaining in atomic form, along with CO, in the gas 
phase. 

H 2 0 and 0 2 axe also not significant coolants for quiescent molecular gas. In the 
case of H2O, a number of known chemical processes can locally elevate its abundance 
in regions with enhanced temperatures, such as warm regions surrounding young stars 
or in hot shocked gas. Thus, water can be a locally important coolant. The new 
information provided by SWAS, when combined with recent results from the Infrared 
Space Observatory, also provide several hard observational constraints for theoretical 
models of the chemistry in molecular clouds and we discuss various models that satisfy 
these conditions. 
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i Subject headings: ISM: molecules — abundances — evolution; stars: formation; molecular 
processes 


1. Introduction 

Since the recognition that reactions between ions and molecules play a major role in 
molecular synthesis (Herbst & Klemperer 1973; Watson 1973; Dalgarno & Black 1976), it has 
long been proposed that H 2 O and O 2 are important species in the dense interstellar medium. 
Indeed steady-state gas-phase chemical theory predicts that H 2 0 and 0 2 , along with O and CO, 
constitute the gas-phase reservoirs of elemental oxygen in molecular gas that is well shielded from 
ultraviolet radiation (e.g., Millar 1990). As a consequence H 2 0 and 0 2 have been suggested as 
important coolants for dense gas (Goldsmith & Langer 1978; Hollenbach 1988; Neufeld, Lepp, & 
Melnick 1995). 

Atmospheric obscuration has impeded direct observation of these two molecules and our 
overall understanding of interstellar oxygen chemistry has remained uncertain. Ground-based 
observations of isotopologs, such as H 2 18 0 and 0 18 0, have provided some information on their 
relative importance (e.g., Zmuidzinas et al. 1995; Marechal et al. 1997) and the Infrared Space 
Observatory has recently provided direct observations of H 2 0 in warm gas (cf. van Dishoeck et al. 
1999). However, despite these efforts, over the past decade laboratory determinations of reaction 
rates and crucial branching ratios has reached the point that our knowledge of the chemical 
pathways theorized to create 0 2 and H 2 0 in cold (T < 50 K) gas has, surprisingly, surpassed our 
observational capabilities. 

The Sub millime ter Wave Astronomy Satellite (SWAS) has opened a new window on 
interstellar oxygen chemistry by conducting pointed observations in low energy transitions of 
ortho-H 2 0 and 0 2 toward numerous dense molecular cores. Some results from these observations 
have been presented in this issue where chemical abundances (and upper limits) relative to H 2 are 
derived. In this Letter we examine the consequences of these results on the chemistry of molecular 
clouds, on the potential of H 2 0 and 0 2 as coolants, and on the oxygen budget in the interstellar 
medium. In our analysis we assume the ortho/para ratio of water is 3:1. 


2. Primary Oxygen Reservoirs 

The total oxygen budget is constrained by observations of diffuse clouds where the amount of 
oxygen not locked inside the refractory cores of dust grains is ~ 6.4 x 10 4 relative to H 2 (Meyer, 
Jura, & Cardelli 1998). About ~ 15% is in gas phase CO (e.g., Sutton et al. 1984) and SWAS 
observations have demonstrated that in most regions, gas-phase water and molecular oxygen 
account for 1% of the available oxygen (Snell et al. 2000a, b,c; Goldsmith et al. 2000). Thus 
gaseous H 2 0 and 0 2 are not primary reservoirs of atomic oxygen. 
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In the solid state, a large portion of the oxygen budget is frozen on the surfaces of cold 
dust grains. Towards NGC7538 (IRS9) ~40% is in the form of H2O, CO, CH3OH, and CO2 ice 
(Whittet et al. 1996). 1 In the cold absorbing gas towards the background star Elias 16 molecular 
ices account for nearly 65% (Schutte 1999). Thus in star forming cores (NGC7538, OMC-1, 
Ceph A, etc.) gas and solid phase molecules account for ~55% of the available oxygen (see also 
van Dishoeck et al. 1999). The remainder could be in the form of atomic oxygen, which has 
been found to be abundant towards a few sources (see Caux et al. 1999 and references therein). 
Another possibility is that there is higher depletion of elemental oxygen into refractory dust cores 
in denser regions as found for other elements (see Snow et al. 1998; Meyer et al. 1998). 


3. Cooling 

There have been several theoretical studies of the cooling from molecular emission (Goldsmith 
& Langer 1978; Hollenbach 1988; Neufeld, Lepp & Melnick 1995; hereafter GL78, H88, and 
NLM95), each focusing on the gas-phase chemistry and resulting cooling for specified conditions. 
In the case of 0 2 , GL78 assume an abundance of x(0 2 ) = 10 -5 and find that 0 2 contributes 
~ 30 - 40% of the total cooling for n Hi < 10 5 cm" 3 and T sas = 10-20 K. NLM95 used a higher 
0 2 abundance of ~ 9 x 10 -5 and found that 0 2 accounts for ~ 20% of the cooling for low density 
gas. Observations of O2 (e.g. , Goldsmith et al. 2000) find x(02) <10 m regions of high and 
moderate density (n# 2 > 10 3 cm -3 ). Since 0 2 emission is optically thin the cooling rate will scale 
with abundance and therefore the contribution of molecular oxygen to the overall cooling must be 
< 1 %. 

Both GL78 and NLM95 assume a water abundance of ~ 10 -6 and for clouds with T sas < 30 
K H2O is predicted to account for rv 10% of the total cooling, provided n Hi > 10 5 cm 3 . In regions 
with low UV intensities, and for T 90 5 > 400 K, a series of neutral-neutral reactions with activation 
barriers will convert all of the atomic oxygen not locked in CO into H2O. Under these conditions 
the abundance of water can be as high as ~ 3 x 10 4 and water is a major coolant for molecular 
gas (H88, NLM95). 

SWAS observations of water in extended gas in Orion and M17 (and other sources) show that 
the abundances in the cold, T < 50 K, quiescent gas are only ~ 10 3 (Snell et al. 2000a, b,c). 
With this abundance, and using the cooling functions provided by NLM95, the total cooling rate 
for o-H 2 0 is a factor of 25 below the rate for x(H 2 0) = 10 -6 . Thus, globally water emission 
is not an important coolant for molecular clouds. However, locally the cooling of water can be 
quite important. The detection of o— H2O emission in a beam size of 3f3 x 4.5 by SWAS does 
not preclude the existence of regions smaller than the SWAS beam with high water abundance. 


‘Assuming the total column density is N(H 2 ) ~ 5.0 x 10 22 cm -3 from the 13 CO observations of Mitchell et al. 
(1990), x(CO) = 2.0 xlO -4 (Lacy et al. 1994), and 12 C/ 13 C = 60. 
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Indeed, ISO observations have shown that water is found in high abundance both in hot shocked 
gas and in warm regions surrounding young stars (see van Dishoeck et al. 1999). For those 
conditions H2O emission is an important, if not the dominant coolant (Saraceno et al. 1999; 
Melnick et al. 2000). 


4. Observational and Theoretical Constraints on Chemical Models 

In warm gas SWAS and ISO have demonstrated that chemical theory is in agreement with 
observations (van Dishoeck et al. 1999; Melnick et al. 2000; Neufeld et al. 2000). For cold 
quiescent gas we view SWAS observations and recent results from ISO as having set several hard 
observational constraints for chemical theory. We suggest that successful chemical models satisfy 
the following conditions: 

1. Produce an abundance of molecular oxygen that is less than 10 -6 (Goldsmith et al. 2000). 

2. Be able to account for the low water abundance (10 -9 < x(H20) < few xlO 8 ; Snell et 
al. 2000a, b,c; Ashby et al. 2000) found in star forming regions with the following physical 
properties (nff 2 > 10 5 cm -3 ; T ffas = 20 — 50 K; T dust £ 50 K; Ay > 10 mag; Go = 1 — 10 5 ) 2 
and for star- less cores (e.g., TMC-1) produce an water abundance < 7 x 10 8 (for n# 2 — lb 
cm" 3 ; T gas = 10 K; T duJt < 15 K; Ay > 10 mag; G 0 = 1). 

3. Be capable of simultaneously reproducing the high water abundances (x(H20) ~ 10 6 ; Neufeld 
et al. 1999) found in the low density gas along the line of sight towards SgrB2 (n# 2 <10 4 cm -3 ; 
T 9 „ s ~ 10 - 30 K). 

4. Account for the co-existence in the gas phase of complex carbon-bearing species (e.g., HC3N, 
HCN, and CH 3 C 2 H) with that of simple and complex oxygen-bearing molecules, such as SO 
and CH3OH (Bergin et al. 1997). 

5. Account for a water ice abundance of x(H20), ce ~ 1 x 10 -4 (Tielens et al. 1991; Schutte 
1999). The appropriate physical conditions are not highly constrained, but are between n ff2 
= 10 3 — 10 6 cm -3 and T 9as ~ 10 — 50 K. 

6. Be consistent with other species (e.g., NH3) observed to have non-negligible concentrations in 
molecular gas and have similar sensitivities to the radiation field. For example, H2O, O2, and 
NH 3 have similar sensitivities to both internal and external radiation fields (van Dishoeck 1988; 
Sternberg & Dalgarno 1996). Therefore the low abundances found by SWAS cannot be the 
result of enhanced UV penetration. 


2 Go = enhancement factor above normal interstellar UV radiation field 
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5. Oxygen Chemistry in Molecular Clouds 
5.1. Time Dependent Chemistry 

The oxygen chemistry begins with the cosmic ray ionization of H 2 which quickly leads to the 
formation of the trihydrogen ion (H^ - ) . H J will react with O and a sequence of rapid reactions 
leads to the formation of the H 3 0 + ion. The next step involves the dissociative recombination of 
H3O+ with electrons: 


H 3 0 + +e~ 


H + H 2 0 fi 

OH + H 2 f 2 

OH + 2H f 3 

O + H + H 2 f 4 


(Rl) 


where /1-4 are the branching ratios. 3 The formation of 0 2 follows through a reaction between O 
and OH. The destruction of both H 2 0 and 0 2 is dominated by the cosmic-ray induced photo-field 
and reactions with H + , He + , C + , and H^ (and C in the case of 0 2 ). 

Fig. 1 presents this current view of pure gas phase time-dependent chemistry for typical 
conditions. The initial conditions simulate diffuse cloud abundances and are taken from Millar et 
al. (1997). This model is somewhat limited in application in that it does not include dynamics and 
the initial conditions are an approximation. Nonetheless it is illustrative of the relative importance 
of the various species and the timescales involved under the assumption of constant density. 

The predicted water abundance changes sharply depending on the assumed H 3 0 + branching 
ratio; for < 0.05 the water is mainly produced via proton-transfer reactions between H 3 0 + and 
other molecular species (e.g., HCN, NH 3 ). Examining the constraints listed in § 4, the model 
is consistent with SWAS 0 2 observations, provided the chemistry for all sources surveyed has 
evolved for t < 10 s yr. The model is also consistent with the high water abundance found along 
the line of sight towards SgrB2, provided fi > 0.3. If, for consistency with the SgrB2 observations, 
we adopt fi = 0.33, then the water abundances found in star forming regions place even greater 
strains on the theory, as the observations require ages < 1000 yr. 4 


3 Vejby-Christensen et al. ( 1997 ) measured /, = 0 . 33 , h = 0 . 18 , f 3 = 0 . 48 , and ft = 0 . 01 , while Williams et al. 
( 1996 ) find fi = 0 . 05 , / 2 = 0 . 36 , / 3 = 0 . 29 , ft = 0 . 30 . These two laboratory measurements differ sharply in the 
amount of water produced. 

4 The ages of molecular clouds and cloud cores are highly uncertain. Lower limits can be set using the free-fall 
timescale, which is ~ 10 5 yr for typical densities (n h 2 = 10 5 cm" 3 ). The mean age of embedded protostars sets 

another lower limit. In Orion the mean age is ~1 Myr (Hillenbrand & Carpenter 2000 ). 
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5.2. Alternate Solutions 

One possible method to lower the water and molecular oxygen abundance is a change in the 
elemental C/O ratio above the solar value of 0.4. If the C/O ratio is raised to higher values (e.g., 
by depleting the elemental oxygen relative to that of carbon) then the amount of atomic oxygen 
available to form water and molecular oxygen is limited and relatively more carbon atoms and ions 
exist to destroy H 2 0 and 0 2 . Additional effects have been discussed in the literature (Terzieva & 
Herbst 1998 and references therein). 

In Fig. 2 we show points defining the plane of 0 2 and H 2 0 abundances as a function of time 
and the assumed elemental C/O ratio. Clearly changing the C/O ratio has a dramatic effect on 
the abundance of both molecules. Consistency with SWAS observations of star forming cores 
(constraints // 1 and $2) requires C/O > 0.9 and has the desired effect of allowing for H 2 0 and 
0 2 abundances to agree at later, more relevant, timescales (t > 10 5 yr). Changing the elemental 
C /O ratio above solar could be the result of variations in the interstellar medium due to different 
formation histories, supernovae or planetary nebula activity, and/or a direct by product of the 
gas-grain interaction. 

In Fig. 3 we show a more complex gas-grain chemical model that demonstrates how the C/O 
ratio could be altered using a chemical mechanism outlined in d’Hendecourt, Allamandola, & 
Greenberg (1985). The first panel shows a gas grain model with no surface chemistry. In this 
model oxygen atoms deplete on grains, remain inert, and their abundance is maintained in the gas 
phase by thermal evaporation. The gas phase oxygen enables a chemistry similar to that shown in 
Fig. 1. In Fig. 3b we show a gas-grain model including two surface reactions (O + gr — >• H 2 0(gr) 
and C + gr -> CH 4 (gr)). Here most oxygen is trapped on grains in the form of water ice, which 
will only evaporate if T dust > 90 K. In contrast, volatile CH 4 is released back into the gas phase 
for Td„ 5 t > 23 K and the carbon chemistry is enhanced (see also Brown and Charnley 1991; Nejad 
& Williams 1992). 

This model, with simple and reasonable assumptions for the chemistry reproduces constraints 
(#1) and (#2) with low abundances of H 2 0 and 0 2 at physically relevant timescales, and also 
enhances the abundances of complex species as exemplified by HC3N (#4). Because the depletion 
timescales are longer at lower densities, the inclusion of grains into the chemistry also naturally 
allows for a difference between the high water abundances found along the SgrB2 line of sight and 
the low water abundances in star forming regions. 

However, there are two potential problems with this model. First, the sulfur chemistry, 
especially the production of SO (which is observed to be reasonably abundant), is effectively halted 
by the gas-grain interaction (Bergin & Langer 1997; Ruffle et al. 1999). Second, observations of 
atomic oxygen in absorption towards similar regions find high oxygen column densities (e.g., Caux 
et al. 1999). If the oxygen absorption features are probing the entire line of sight (and not just 
cloud halos), then this would be in conflict with the depletion of oxygen in the form of water ice 
(Fig. 3b). 
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Both of these problems could be alleviated by the presence of significant unresolved chemical 
structure along the line of sight with pockets of undepleted molecular chemistry, such as might 
exist at the low density cloud edges or in the vicinity of embedded sources, coexisting (within the 
SWAS beam) with regions of high depletion. The spatial resolution of the SWAS instrument is 
insufficient to resolve such structure. 

Other possibilities that must be considered axe that shocks can alter the C/O ratio by 
creating water ice mantles, supplementing the products of surface chemistry (Bergin et al. 1998), 
bistability in chemical models (Le Bourlot et al. 1993; see discussion in Goldsmith et al. 2000), 
and/or that the apparent chemical youthfulness of dense cores could be maintained by turbulence 
or cycling of material between low density to high density phases (Chieze k Pineau des Forets 
1989; Charnley et al. 1988; Xie, Allen, k Langer 1995). Indeed in these models the 0 2 and H 2 0 
abundances can be lowered several orders of magnitude below steady state values as required by 
SWAS observations (Chieze k Pineau des Forets 1989; Xie, Allen, k Langer 1995). 

We acknowledge valuable suggestions given by an anonymous referee. This work was 
supported by NASA’s SWAS Grant NAS5-30702 and by grant NAG5-3542 from NASA’s Long 
Term Space Astrophysics Program. 
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Fig. 1. — Time evolution of chemical abundances relative to H 2 for a solar C/O ratio of 0.4. The 
evolution of the H 2 O abundance is shown for different values of the H 30 + branching ratio (discussed 
in § 5). The assumed physical conditions are: n# 2 = 10 5 cm"" 3 , T gas = 30 K, and Av = 20 m . The 
chemical model used here has been described previously (e.g. Bergin & Langer 1997) and is similar 
to other current gas phase models (Lee et al. 1996). 
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Fig 2 . — Plot of O 2 versus H 2 O abundances shown as a function of time and assumed elemental 
C /O ratio. The different point types delineate the various C/O ratios (a value of 0.4 is the solar 
ratio). The solid grey lines sure isochrones of constant time. Models are run using n h 2 ~ 10 5 cm » 
T 9as = 30 K, and Ay = 20 m . The shaded surea denotes the range of abundances implied by SWAS 

observations. 
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Fig. 3 — Abundances of CO, 0 2 , H 2 0, and HC 3 N as a function of time for two models of the 
gas-grain chemistry, (a) Model including depletion and desorption onto grain surfaces, but no 
surface chemistry, (b) Model including depletion and desorption onto grains but with simple 
surface chemistry. Both models assume n# 2 = 10 5 cm -3 , T 9as = Td ust = 30 K, and Ay = 20 m . 
The second model does not include any modifications to account for the discrete nature of grains 
(Caselli, Hasegawa, & Herbst 1998), but is nonetheless illustrative of the effects of the process. 




